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The Schrtidinger functional is used to define a renormalised coupling for pure SU(4) Yang-Mills 
theory, with Wilson action and suitably selected boundary conditions on the link field. The cou- 
pling, which runs with the size of the lattice, is then determined by a finite-size scaling technique 
through a large range of momenta, thereby allowing a connection to be made between the high 
energy regime and the low energy non-perturbative regime, where contact is made with the phys- 
ical scale of the theory. Using data from previous SU(2) and SU(3) simulations obtained with the 
same technique, the running of the 't Hooft coupling defined through the Schrodinger functional 
is studied, and we check whether the large-rV expectation that g 2 N is a universal function of the 
energy scale E holds down to energies of the order of the string tension. Finally, we determine 
A-SF I \[g as a function of N at leading order in 1 /N 2 . 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



Running coupling in pure SU(4) Yang-Mills theory 



Gregory Moraitis 



1. Introduction 

In perturbative calculations in QCD, dimensional transmutation is accounted for by a scheme- 
dependent quantity with mass dimension one called the A-parameter. Since QCD is a one-scale 
theory, A must be related to the non-perturbative scale of the theory, which can then be used to 
express the perturbative running of the coupling. In this respect, the lattice has proved to be a 
valuable tool, and the running coupling has already been computed for pure SU(2) and SU(3) 
Yang-Mills theories [1, 2] using a scheme referred to as the Schrodinger functional [3]. Here, we 
apply this method to the SU(4) theory. 

In the Schrodinger functional scheme, the renormalised coupling runs with the physical size 
of the system L. To cover a wide range of energies, an iterative step-scaling method is used [4]. 
A step function o{s,g 2 {L)) = g 2 (sL) is introduced and determined pointwise for systems of size L 
and sL in every direction. To step up the scale, the bare coupling is tuned from go to g' , such that 
g 2 (L)\ g i = g 2 (sL)\ g0 . This renormalisation group transformation can be performed iteratively and 
effectively allows us to accommodate larger system sizes on fewer lattice sites. Of course, once on 
a lattice, the system is affected by the finite lattice spacing a so, in practice, each physical box size 
is simulated for different lattice spacings and an extrapolation to the continuum is performed. With 
this method, we successfully cover a range from 7GeV down to 0.44GeV. 

At the lower end of this range, contact can be made with the non-perturbative observables of 
the theory. In particular, by using the data in [5] and the interpolating formulae in [6], we are able 
to express the coupling in physical units. Once the scale is fixed, we can obtain a value for A$f in 
the Schrodinger functional scheme, and compare the lattice data with the perturbative evolution of 
the beta function to two-loops. 

2. Method 

2.1 Preliminaries and notation 

The basic ideas and method are taken over directly from [3], so only the basic notation is 
outlined below. 

The Schrodinger functional 3? is defined on the lattice with the standard Wilson action, 

3f[C,C'} = jD[U}e- s ^ u \ (2.1) 

S[£/] = ~£Tr(l ~U(p)), (2.2) 
Sop 

where U (p) denotes the plaquette, and the sum must be explicitly taken over both orientations. 
The boundary links are required to satisfy inhomogeneous Dirichlet boundary conditions for 

it = 1,2,3, 

W(x,k)\ xQ=0 = exp(aC k (x)), W(x,k)\ x x, =L = exp(aC[(x)), (2.3) 

where Q and C' k are spatial boundary fields which need to be chosen (cf. §2.2). 

The fixed boundary conditions induce a colour background field B into the system so that the 
effective action can be written as an asymptotic series 

r[B] = gQ 2 r [b] + n [b] + g 2 r 2 [b] + ■ ■ • . (2.4) 
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By introducing a dependence of the boundary links (and thus of the background field) on a real 
dimensionless parameter rj , we can then define a renormalised coupling as 



(2.5) 



for a particular choice of v\ . On the lattice, this quantity can be measured via Monte Carlo simula- 
tion, by calculating the expectation value of the observable 



= - JEE + E 'iW + +^'( x 



))■' 



Ei (x) = Tr [cW(x, I) U (x + at, 0) U (x + aO, I) 1 U (x, 0) f ] x0=Q , 



(2.6) 



where c is a matrix appearing when differentiating the boundary field dependence of the Wilson 
action (cf. §2.2), and a similar expression holds for E'(x). 

2.2 The fundamental domain in the SU(4) theory 

It is desirable to choose boundary fields which minimise the effect of the finite lattice spacing. 
It was shown in [3] that, for N colours, the optimal choice are constant Abelian fields 
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(2.7) 



Stability considerations of the background field constrain the angles 



N 

I 

i=\ 



0/ = 0, 0i < 02 < - < 0jv, |0i - 0/1 < 2^ 



(2.8) 



and similarly for 0' (from here on we drop the suffix k on the angles and use the same choice for 
k = 1,2,3). A set of angles satisfying these constraints is said to be in the fundamental domain. 

Specialising to SU(4), the fundamental domain can be described symmetrically by defining 
a one-to-one map between the set of angles (01,02,03,04) and a point V in a certain bounded 
three-dimensional region, 



V= T(0i-ei + 02-e2 + 03-e 3 + 04-e 4 ) 



0; = V-e/, 



(2-9) 



where e, are the weights of the Lie algebra of SU(4) in the fundamental representation, normalised 
as e, ■ ej = I (4<5,-y — 1 ) . With this normalisation, the vertices v, of the fundamental domain are 



3k 



E e 7' 



i= 1,2,3,4, 



(2-10) 



7=1 



describing a skewed tetrahedron (figure 1). Having identified the fundamental domain, one must 
then select the set of angles for each boundary, Q and C' k , each set corresponding to a point inside 
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Figure 1: The fundamental domain and its two planes of symmetry. The points marked x correspond to the 
angles (2.1 1). 

the tetrahedron. These points will be members of a one-parameter family of angles, parameterised 
by TJ, through which the renormalised coupling is defined by (2.5). In principle this choice is 
arbitrary and of no conceptual significance, however the signal-to-noise ratio of the Monte Carlo 
evaluation is highest when (i) the points are well away from the domain's edge, (ii) they are as far 
as possible from each other, and (iii) the two boundaries are on an equal footing. Geometrically, 
this corresponds to choosing two points related by a symmetry of the fundamental domain. 

Keeping the above points in mind, we choose the two points to be related by the symmetry 
reflecting about the plane through V2 and V4, and make the particular choice 

01 = -51? - ?wV2 ^ = i|]-i(2 + V2)ff 

02 = 4r]-i(2-V2> 02 = h-?( 4 -V2)* 

03 = iT] + i(2-\/2> $ = -iT] + A(4-x/2> 

fa = fa + faV2 $' A = -\i} + \(2 + y/2)%. (2.11) 

With this choice, we set T] = to obtain 

r'M = -—r^[^\- (2.12) 



a 1 \2L 2 , 

The matrix c defined in (2.6) is then c = diag(— A, — A, A). 
3. Results 

The results of our simulations are summarised in table 1. The data is grouped into four blocks, 
each corresponding to a fixed value of g 2 (L). We use a step of s = 2 throughout, and calculate the 
renormalised coupling g 2 (2L) for different lattice spacings, extrapolating to the continuum. The 
eiTors quoted are statistical only, and at this stage all results should be regarded as preliminary. 

To convert this dimensionless scale (L/a) to a more physical one, we need to set the lattice 
spacing at a fixed value of /3. The quantity a^fo is readily available in the range 10.55 < /3 < 
11.40 [5], so we can deduce the physical lattice spacing at the highest couplings probed, ie. at 
g 2 (L) = 2.46 (table 2). Using the interpolating formula in [6], we obtain a^fo = 0.16069 at j3 = 
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J8 L/a a/L g 2 (L) g 2 (2L) 


Extrapolation 


15.126 6 0.167 1.0222(6) 1.2247(10) 
15.626 8 0.125 1.0223(7) 1.2162(18) 
16.000 10 0.100 1.0223(5) 1.2104(21) 


1.1891(37) 


14.137 6 0.167 1.1893(3) 1.4832(13) 
14.634 8 0.125 1.1894(4) 1.4700(19) 
15.007 10 0.100 1.1890(8) 1.4634(25) 


1.4332(47) 


13.142 6 0.167 1.4328(3) 1.9042(14) 
13.631 8 0.125 1.4332(4) 1.8797(24) 
14.000 10 0.100 1.4330(6) 1.8667(35) 


1.8099(62) 


12.190 6 0.167 1.8098(7) 2.7198(38) 
12.668 8 0.125 1.8102(6) 2.6548(49) 


2.4601(300) 


Table 1: Pairs of running couplings for fixed values of j3 = 8 /go- 



L/a 


P 


6 


11.329 


7 


11.571 



Table 2: Bare coupling vs lattice sites at fixed g 2 (L) — 2.46. 

11.329 and a^fa = 0.13523 at /3 = 11.571. These two values are used individually to compute the 
scale of the theory and the results match to better than 2%. Using ^fc = 420 MeV, the maximum 
(minimum) box size is 0.45fm (0.028fm) corresponding to an energy range from 0.44GeV to 7GeV. 
From asymptotic freedom, we expect that perturbation theory holds accurately at the upper end of 
the scale, so we impose that g 2 (7GeV) = 1.0222 and use the two-loop solution to the RG equation 
for the coupling g 2 at an energy scale E 

Asf VA) Pog 2 (E)J 
to obtain a value for the A-parameter for N = 4 colours, 

A SF (N = 4) = 101 ±1 MeV. (3.2) 

The error quoted only reflects the scale uncertainty which comes from fixing the lattice spacing at 
the two values of j8 , though this is expected to be the largest source of error. 

Figure 2 shows the lattice data plotted alongside the perturbative evolution predicted by the 
beta function. The data is thus accurately described by two-loop perturbation theory down to energy 
scale of the order of ^fa . 

Scaling with N 

It is expected that for large -A^ the running of the t'Hooft coupling g 2 N is independent of N, 
so that a universal curve describes the running for any N. Previous lattice simulations have hinted 
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Figure 2: Results of lattice simulations for the running coupling in the SU(4) theory, together with the one 
and two loop perturbative predictions. 



that this holds also at energy scales where one would expect large -N perturbation theory to break 
down [7]. We are now in a position to test this using the results herein, together with those of SU(2) 
and SU(3)in [1,2]. 

As for the case of SU(4), we fix the physical scale of the SU(2) and SU(3) data using updated 
values of the string tension [5,6], and superimpose the results for /V = 2, 3,4 (figure 3). The points 
fall on a single curve within the errors, suggesting that the prediction of universality is accurate 
even at N = 2. Again, the errors shown correspond to the scale uncertainty introduced when fitting 
to the string tension and we expect this to be the biggest source of error, especially in the case of 
SU(2) where the range of /3 for which we have accurate values of the string tension does not overlap 
with the range of j3 covered in the Schrodinger functional simulation [1]. The corresponding A$f 
were found to be 



A SF (N = 2) = 110±10MeV, A SF (N = 3) = 103±3MeV. 



(3.3) 



From general theoretical considerations, we expect that the leading order corrections to the univer- 
sal behaviour are of order l/N 2 . In fact, our data are consistent with this expectation. Though a 
fuller analysis of the errors should be carried out at this stage, we tentatively find that, in the range 
of N considered, A$f(N) is well approximated by 



A SF (N) 



■ 0.23 + 



0.12 



(3.4) 



4. Conclusions 

We have generalised the Schrodinger functional to SU(4), and used it to produce data for the 
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Figure 3: Running of the t'Hooft coupling g 2 N for N = 2,3,4. 



running coupling from 7GeV down to the scale of ^/a, effectively linking the perturbative to the 
non-perturbative regime. The result is similar to the previous SU(2) and SU(3) calculations already 
performed in that the data is well approximated by perturbation theory down to energies of the order 
of the string tension. Of course, this should not be taken to mean that perturbation theory correctly 
accounts for all phenomena down to this energy, as the coupling defined through the Schrodinger 
functional may simply be an exceptional case. 

Similarly, the t'Hooft coupling g 2 N is seen to be a universal function of E, only weakly de- 
pendent on N even down to N = 2, and we have extracted the N dependence of A$f to leading 
order in l/N 2 . The large -N expectation of universality thus holds true down to energies of order 
yfo, though we are unable to make any certain claims as to whether universality really extends to 
the non-perturbative level, since perturbation theory itself proved adequate in the ranges of energy 
studied. 
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